Multiple bound states in scissor-shaped waveguides 
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We study bound states of the two-dimensional Helmholtz equations with Dirichlet boundary condi- 
tions in an open geometry given by two straight leads of the same width which cross at an angle 0. 
Such a four-terminal junction with a tunable 6 can realized experimentally if a right-angle structure 
is filled by a ferrite. It is known that for 6 — 90° there is one proper bound state and one eigenvalue 
embedded in the continuum. We show that the number of eigenvalues becomes larger with increas- 
ing asymmetry and the bound-state energies are increasing as functions of 9 in the interval (0, 90°). 
Moreover, states which are sufficiently strongly bent exist in pairs with a small energy difference 
and opposite parities. Finally, we discuss how with increasing the bound states transform into the 
quasi-bound states with a complex wave vector. 

PACS numbers: 42.25.Bs,03.65.Ge, 73.40.Lq,03.75.Be 



I. INTRODUCTION 

The question of a possible existence of modes trapped 
in open two-dimensional systems has been a classic in the 
theory of waveguides; trapped modes due to particular 
boundary conditions were studied already half a century 
ago However, only much later it was realized that 
the introduction of bends and crossings into waveguides 
gives rise generally to confined states, or bound states 
which exist below the cutoff frequency for the waveguide 
1, 1, |, 1, 0, g 1, 0, H 0. The 

existence of such states has both theoretical significance 
and implications for possible applications. They have 
been subsequently discussed in many papers, in addition 
to those mentioned above we refer to the monograph |l7j 
and the bibliography there. 

In this paper we consider a system of two straight 
waveguides of the same width d which cross at a nonzero 
angle 8. The right-angle case was of one of the first ex- 
amples where the binding was studied. Schult, Raven- 
hall, and Wyld showed the existence of two bound 
states. One of them is a true bound state at the energy 
0.66(7r/<i) 2 in the natural units, while the other one at 
3.72(7r/ii) 2 is embedded into the continuum and does not 
decay due to the symmetry. The latter corresponds to 
the single bound state in an L-shaped tube of width d/2 
Q . Our aim is to show how the spectrum of such a junc- 
tion, which we will call for the sake of brevity "scissors" 
in the following, changes as the angle 8 varies over the 
interval (0,90°). 

We will show that as we go further from the cross 
symmetry if the right-angle structure, new bound states 
emerge from the continuum. In strongly skewed junction 
corresponding to a small 8 there are many of them. The 
mechanism responsible for their existence is the same as 



for the bound states in sharply broken tubes studied the- 
oretically and experimentally in |jTfj[ |L2| , namely a long 
part of the junction where there transverse contribution 
to the energy is substantially lower than (n/d) 2 . In the 
present case, however, the system has a mirror symmetry 
with respect to the axis of the complement angle 180° — 8 
and the bound states exist in pairs corresponding to dif- 
ferent parity. We will show that as the angle 8 diminishes 
and the states become strongly bound, the energy gap 
between the even and odd member of the pair vanishes 
exponentially fast. We also study the behavior around 
the critical values 8 C where the bound states emerge from 
the continuum. Our numerical analysis shows that the 
binding energy of the weakly-coupled states behaves as 
correction w ir 2 — r y(8 c — 8) 2 



II. FERRITE FILLED MICROWAVE 
WAVEGUIDES AS A WAY TO VARY THE 
ANGLE OF THE SCISSORS 

Before analyzing the scissor spectrum, let us discuss 
how such a structure can be realized experimentally as a 
microwave device. It is no problem, of course, to build 
crossed waveguides in the ways explained in |l7| ]. How- 
ever, in such a setting it is not easy to vary the geom- 
etry continuously. Our point here that this goal can be 
achieved with a structure of a fixed angle if the latter is 
filled by a ferrite with an axial magnetic anisotropy and 
an external magnetic field is applied. We will show that 
this leads to an effective angle controlled by the field 
strength, following the idea which was first applied to 
the equivalence between a ferrite-filled squared resonator 
with an external magnetic field and a field-free rhombic 
polygon E§). 
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To explain the mechanism of this equivalence we begin 
with the Maxwell equations which in the presence of a 
material have the form 

V-E = V-B = 0, 

V x E = -ikB, V x H = ikE, 

B = pM, (1) 

where E is the electric field, H is the magnetic field, B 
is the magnetic induction, k = lo/c and u> is an eigen- 
frequency with the wave vector k. We suppose that the 
material has a magnetic anisotropy corresponding to an 
anisotropic permeability p, = 1 + 4-7TX with |n| 



(Xxx Xxy 
Xxy Xyy 
0, 



(2) 



where 



and Eq(3) 



gfhMo 
fiiOa — u> 2 ' 



Xxy Xyx 



Xyy 



^1^2 — U 2 



iujgCliMo 



2 ' 



(3) 



Here g is the magnetomechanical factor, Mo is the mag- 
netization of the material, 

, r ( M H^ gK a , \ 
fix = gM I , J + ^ cos 2 * 



(i) 



tt 2 = gM 



M H - 

M 2 



Mn 



gKa 

Mn 



cos 2* 



(4) 



and K a characterizes the anisotropy type: it is an easy 
plane anisotropy for K a > and an easy axis anisotropy 
for K a < 0. In what follows we suppose that the material 
has an easy plane magnetic anisotropy, K a > 0, in which 
case the intrinsic magnetic field is equal to 

H = Ho-47tM 0z . 

On the relations @) ^> is the angle between the anisotropy 
axis N and the magnetization Mo. We choose the latter 
to coincide with the z axis along the magnetization, while 
the x axis lies in the plane spanned by the vectors N and 
M . 

In the simplest case of an easy plane magnetic material 



we have M _L N and M 
obtain from (Q) 



with * = tt/2, so we 



Q x = g(H a - iirMo), 

r» 2 = g(H - iirMo) + gK a M a . 



(5) 



This is shown in Fig. [I] where the z axis is perpendicular 
to the plane of the waveguide. We seek a two-dimensional 
solution of the Maxwell equations ([!]) shown at this figure 



FIG. 1: Schematical view of the cross-bar resonator (scissors 
with 9 — 90°) filled with ferrite where M is the magnetization 
of ferrite and N is the anisotropy field. 



in the form E(x,y) = E(x,y)e z . The fields B,H lay in 
the plane x, y and depend on x,y too. Then the first 
equation is satisfied, while the third one gives 



— ikB r = 



dE z 
dy ' 



-ikB,, 



dE z 
dx 



(6) 



and finally, the fourth Maxwell equation can be rewritten 
as 



ikE, = 



dx 



dH x 
dy 



{<) 



Using the explicit form of the permeability given by 
we get 

(\^xx \^xy 
(J>yx Hyy 
1 

Combining this with Eq. (^J) we obtain 




H, 



j_ ( f^yy f^xy 



i dE z 
k dy 
-i dE z 
k dx 



(8) 



where we have denoted 

■D l^xxf-Lyy f^xyf^yx ■ 

Substituting Eq. (|) into Eq. @ we obtain 



d 2 E z 

dx 2 



d 2 E z 

dy 2 



+ (pxy+»yx)-^£+Dk 2 E z =Q. 

(9) 

The key observation is that the mixed derivatives in the 
last equation can be eliminated by the coordinate trans- 
formation 



Ma; y ~\~ j-^yx 




(10) 
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which allows us to cast Eq. (£)() into the following simple 
form 



V 2 E Z + Xk 2 E z = 0, 



where 



f^XX 



(11) 



(12) 



we have taken into account that /i xy + /j, yx = holds in 
accordance with (||). 

The transformation (|l0|) defines a relation between a 
right-angle cross structure and a skewed one with an an- 
gle defined by A. It is too daring, however, to speak about 
a full equivalence, because it is clear from the formulas 
expressing the elements of that the angle depends on 
the eigenfrequencies involved. Let us ask under which 
conditions this dependence of the geometrical factor ( O ) 
can be suppressed. Substituting into ([l2]) the expressions 
© we get 



A 



ftifi 2 - uj 2 + 4Trgtt 2 Mo 



(13) 



Following to JUI we can simplify this expression if to 
assume that 



gK a /M » max{g£fo,47rgM , uj}. 



(14) 



For typical ferrites K a ~ 10 6 erg/cm and AttAIq ~ 
100 Oc |L9]]. Taking the magnetomechanical factor g ~ 
10 7 sec~ 1 Oe~ 1 we obtain from@ #o < 10 4 Oe, uj < 
10 what would require very wide waveguides of a width 
d ~ 10 cm. However, there are ferrites with K a ~ 
10 s erg/cm 3 which lead to the inequality uj <C 10 13 . 
Hence in this case we are able to use standard waveg- 
uides the width of which is of order 1 cm. Then we can 
simplify the geometrical factor of the waveguide to the 
form 



A = 



Ho - AwMq 



(15) 



This formula gives a remarkable possibility to change the 
angle of the scissors 



9 = 2 arctan A 



(16) 



by means of an external magnetic field applied along the 
magnetization direction. 

Moreover if to apply strong magnetic field gHo 3> uj 
or Hq 10 4 Oe, then it follows from formula (|l3|) that 



H 



H Q + gK a /Mo 



(17) 



On this case also the effective angle of the structure can 
be tuned by variation of the external magnetic field. 



III. BOUND STATES AND RESONANCES 

First we review some general properties of the bound 
states which can be derived by the standard methods of 
dealing with Laplace operator with Dirichlet boundary 
conditions, the Dirichlet-Neumann bracketing 20 and 
eigenvalue variation with respect to the change of the 
domain pi] . The problem has two mirror symmetries 
with respect to the axis of the angle 9 which we call 
scissor axis, and with respect to the axis of the larger 
angle 180° — 9 which we call the second axis. Using the 
mentioned methods together with an eigcnvaluc-in-thc- 
box estimate analogous to that employed in jl0| we find 
that 

(i) every bound state is even w.r.t. the scissor axis, 

(ii) with respect to the second axis the bound states 

can have both parities which are alternating if the 
bound states are arranged according to their ener- 
gies, 

(iii) as 9 becomes smaller new bound states emerge from 
the continuum. The number TV of bound states 
satisfies the inequality N > 2ctt~ 1 (90° / 9) with c = 
(l-2~ 2 / 3 ) 3 / 2 « 0.225. While it is not good around 
9 = 90°, where we know that N = 1 from 0, it is 
asymptotically exact as 9 — > 0, 

(iv) all the bound-state energies are monotonously in- 
creasing functions of 9. 

The angle dependence of the bound-state energies has 
different regimes. In the weak-coupling regime when the 
scissors are closing and just passed the critical angle 9 C at 
which a new bound state appeared, our numerical analy- 
sis shows that the binding energy of the weakly-coupled 
states behaves as w ir 2 — r y(9 c — 9) 2 with some constant 
7 which depends on the particular state. On the other 
hand, strongly bound states corresponding to a small 9 
are in the leading order determined by the one dimen- 
sional potential well given by the lowest transverse eigen- 
value ]2^] . The second axis determining the parity of the 
solution is then deep in the classically forbidden region, 
so we can conclude that 

(v) as 9 becomes smaller the bound states group into 

pairs with opposite parities and the energy gap be- 
tween them is exponentially small as 9 — > 0. 

After these general results let us pass to the numerical 
solution. We use three different methods. The most com- 
mon among them is the boundary integral method p3[ . 
In combination with the above general results, it pro- 
vides a rather complete information about the discrete 
spectrum. 

On the other hand, the boundary integral method tells 
us nothing about the scattering problem in the scissor 
structure. We are interested in particular in the scatter- 
ing resonances associated with quasibound states, which 
are characterized of complex values of energy at which 
the analytically continued resolvent has a pole singular- 
ity. A suitable method to treat this problem is the exte- 
rior complex scaling. The method was suggested in the 
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seminal paper p3 and has developed into an efficient 
computational tool - see[^5| and reference therein. The 
use of exterior complex scaling for waveguide structures 
was first proposed in [^6| , here we employ it in the form 
presented in 27 . Before the proper scaling we pass to 



right-angle scissors by means of the coordinate change 



x — x sm a — y cos a 
i 



y 



(18) 



which takes the Hamiltonian to a unitarily equivalent 
operator acting as 



= I - 



d 2 



d 2 



dx' 2 dy r 



2 cos a 



d 2 



dx'dy' 



(19) 



Now we apply the scaling transformation to the longi- 
tudinal variable in the structure arms which leaves the 
central area unchanged, x = g(X) and y = g(Y), which 
yields the scaled Hamiltonian 



H = -V 



with 



c u (X,Y) c 12 (X,Y) 
c 21 (X,Y) c 22 (X,Y) 



■U(X,Y) (20) 



cu(X,Y)-^ 2 
c 2 i{X, Y) 



1 



9' 2 {X) ' 
cos a 



ci 2 (X,Y) 



9'{X)g'(Y) 
1 



and 



U(X,Y) 



9'{X)g'{Y) ' 
2g'{X)g"'{X)-bg" 2 {X) 



Cl2 {X,Y) = — 



9 ,2 (Y) 



4.9'PO 4 

2g'{Y)g"'(Y)-hg" 2 {Y) , g"(X)g"(Y) 



Ag'{YY Ag' 2 (X)g' 2 (Y) 
The function g{x) can be chosen, e.g., as 



2 cos a 



x if \x\ < xq 

6f{x) if \x\ > x 



with xq larger than the channel halfwidth and the inter- 
polating function f(x) such that f(x) — x for \x\ > 2xq, 
the function g{x) is three times differentiable and the 
inverse map g -1 exists. As long as as the parameter 9 
is real, the above transformation is a simple coordinate 
change which does not change the spectrum. However, if 
9 assumes complex values we observe a different behavior 
in the discrete and continuous part typical in such situa- 
tions |^0) : each branch of the continuous spectrum of the 
operator (|20|) is rotated into the complex plane giving 



! (0,oo)} 



71 = 1 



for d = 1. If 3m 9 > the rotated branches point to the 
lower half plane and reveal parts of other sheets of the 






FIG. 2: Bound state energies for scissors structure as a func- 
tion of the interior angle 9. The complex scaling method data 
are shown by points. The boundary integral method data are 
shown by circles. The asymptotic formulas (^) with cor- 
responding quantum numbers m = 1,2,3 are given by thin 
solid lines. Insets above show blow up of asymptotic behav- 
ior of the bound state energies in the vicinity of bottom of 
propagation band pi 2 . 



Riemann surface of energy and we are able to see the res- 
onance poles as complex eigenvalues of the transformed 
operator; the corresponding eigenfunctions are after the 
transformation decaying at large distances, instead of the 
original growing oscillations typical for Gamow functions. 

Finally, the third method is based on application 
of infitesimally small time-periodic perturbation. The 
bound states with energies below the propagation sub- 
band (Eb < Eq = it 2 ) do not participate in stationary 
transmission. However, it is possible to mix the bound 
state | b > with propagating state \k > via a time-periodic 
perturbation 



V(t) = V Q cos(wt) 



(21) 



provided that the matrix elements of the perturbation 
< b\V\k >^ 0. Such a possibility was demonstrated 
for the four-terminal's Hall junction for electron trans- 
mission effected by a radiation field |l5|, psf . Later the 
mixing of bound states with propagating modes was also 
realized in microwave transmission Here we simi- 

lar to 1 29 use the time-periodic perturbation ( pl| ) as a 
probing one to find the bound state energy by resonant 
features in the transmission probability. 



IV. NUMERICAL RESULTS 

Let us show the results of numerical analysis based on 
the methods described above. First we plot the bound 
state energies as a function of the scissor angle 9. The 
results of complex scale method are presented in Fig. || 
by points. 

The results of boundary integral method are shown in 
Fig. B by circles. For the limit 9 — > the energies of 
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FIG. 4: The probability of transmission through the scissors 
structure as function of E = Xk 2 , eigen value of the Helmholz 
equation (|ll|). The transmission probability limits to zero for 
E — *• Eo — tv 2 . Circles show for which E the transmission 
probability equals the half. 



FIG. 3: The even-even and odd-even bound states for scissors 
structure for 6 = 30°. 



bound states are derived in 
form Eq(22) 



O and have the following 



E nm -» — K + (2« 2 + m 2 /A)8 2 l A + 



(22) 



where the quantum numbers n,m = 1, 2, 3, .... The fac- 
tor 1/4 in ( p2| ) takes into account that the length of an 
inscribed rectangle has a length twicely more than in 
the case of bent waveguide studied in JlJ, Insets 
above the figure show blow up of asymptotic behavior 
of the energies in the vicinity of the edge energy band 
Eq = it 2 . For all energies of bound states the asymp- 
totics are ir 2 — r y(8 c — 8) 2 where 7 is a constant. 

As it was discussed in the Introduction as the angle 
8 diminishes and the states become strongly bound, the 
energy gap between the even and odd member of the pair 
vanishes exponentially fast. In fact, one can see in Fig. 
U that the second bound state energy approaches to the 
first one, the fourth bound state energy approaches to 
the third one, and so on. In Fig. || (a, b) the first even- 
even and the second odd-even bound states are shown. 
One can see that we have typical quantum mechanical 
task of double well potential |3l| in which an energy dis- 
tance between the first and second energy levels becomes 
exponentially small with growth of a potential barrier be- 
tween wells. Fig. H (c,d) demonstrates the next pair of 
the bound states in the scissor structure. 

Usually the bound state transforms into the quasi- 
bound state which displayes as resonant dip or peak for 
transmission through the structure as it was observed, 
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FIG. 5: Evolution of resonant features for the transmission 
through the scissor structure caused by mixing the propaga- 
tion state with the bound state via the time-periodic per- 
turbation as the angle of the scissor 6 increases. The 
evolution exactly follows to a blow up shown in inset of Fig. 
^| for 8 — > 8 C — 0. For the angle exceeding 9 C resonant fea- 
ture is decayed. Dashed vertucal line shows the edge of the 
propagation band Eo = n 2 . 



for example, in pOfl . However as one can see from Fig. ^ 
numerical computation of the transmission through the 
scissor's structure does not show any resonant features 
for 8 > 9 C w 71.5°. Also we used the time-periodic per- 
turbation method to search the quasi boun+ state above 
8 C . The results of computation are in the vicinity of crit- 
ical angle 8 C are shown in Fig. M, 

One can see that for 8 < # c we have resonant features 
caused by mixing of propagating mode with the bound 
state of the time-periodic perturbation. However for 9 > 
8 C these resonant features do not evolve but only decay 
with increasing of the angle 8. The small wiggle around 
the value 9.8704 is an artefact of the computation which 
diminishes with the decrease of Vq. 
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FIG. 6: Distances to the bottom of propagation band at which 
the transmission takes the half (see Fig. ^) versus the angle 
of the scissors waveguide. 



Moreover for 9 approaching 9 C the transmission prob- 
ability T(E) undergoes the following feature as shown 
in Fig. |j. The less a value \6 — 6 C \ the more a slop of 
the transmission versus E in the vicinity of E — 7r 2 . For 
9 -> 9 C the derivative dT/dE -> oo. If to plot E - n 2 
at which the transmission is reaching half as shown by 
circles in Fig. ^ versus the angle of scissor waveguide we 
obtain remarkable curve shown in Fig. ^ We see that a 
value E — it 2 equals zero just at the critical angles. 
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